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The origin of the correlation between surface reactivity and quantum-size effects, observed in 
recent experiments on the oxidation of ultrathin magnesium films, is addressed by means of ab initio 
calculations and model predictions. We show that the decay length in vacuum of the electronic local 
density of states at the Fermi energy exhibits systematic oscillations with film thickness, with local 
maxima induced when a quantum well state at feii = crosses the Fermi energy. The predicted 
changes in the decay length are expected to have a major impact on the electron transfer rate by 
tunneling, which has been proposed to control the initial sticking of O2 in the oxidation process. 
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There is considerable interest in identifying methods 
to tailor the chemical reactivity of surfaces, a crucial fac- 
tor in many technologically relevant surface phenomena, 
including oxidation and catalysis. The recent observa- 
tion of a correlation between quantum-size effects and 
the surface reactivity of ultrathin metal films Q is an 
exciting development in this area, both because of its 
fundamental interest for a quantitative understanding of 
the structure-size dependence of the reactivity and in 
view of the importance of nano-scale and low-dimensional 
structures in modern technology. Experimentally, oscil- 
lations were observed in the oxidation rate of ultrathin 
Mg(0001) films on W(110), as a function of film thick- 
ness, with the largest oxidation rate occurring when a 
quantum-well state was found to cross the Fermi energy, 
in photoemission (PE) spectra taken near normal in- 
cidence. In particular, the changes observed in the initial 
oxidation rate — when most of the film was still metallic, 
were found to be dramatic. The origin of these changes 
has not been established yet. The changes in reactivity 
were suggested to be due to a change in the density of 
states (DOS) at, or near the Fermi energy Ep Q], as the 
oscillations in the reactivity were found to correlate with 
the oscillations in the PE intensity at Ep. However, un- 
like the normal-incidence PE intensity, which measures 
the partial DOS at fc|| = 0, the total DOS of the films 
at Ep is bound to exhibit a monotonic, staircase-like in- 
crease with film thickness 0, and cannot thus simply 
account for the oscillations in the reactivity. 

Here we propose a theoretical interpretation of these 
observations, in which the decay length in the vacuum of 
the electronic local density of states around the Fermi en- 
ergy is identified as the key parameter responsible for the 
changes in the reactivity. On the basis of ab initio calcu- 
lations, inspired by model predictions for quantum well 
states, we show that the decay length, A, of the local den- 
sity of states at the Fermi energy exhibits an oscillatory 
behavior with film thickness, with local maxima present 
when the highest occupied quantum-well state at fc|| = 



is closest to the Fermi energy. The predicted changes in 
the decay length are significant. They are expected to 
have a direct, exponential impact on the electron trans- 
fer rate by tunneling, which has been proposed to control 
the initial sticking of O2 on the metal surface The 
changes in A should be observable by scanning tunneling 
microscopy (STM), near terrace steps, on surfaces char- 
acterized by regions of different film thicknesses . 

Model predictions for A. — Within an independent elec- 
tron description and considering, for simplicity, a square- 
well potential along the z direction, normal to the thin 
film, with a constant depth V and a variable width L 
(film thickness), the electronic states are solutions of a 
separable problem in the z and (x, y) variables, and read: 

^ M (x,y,z)^ X n(z)-e i ^+ k ^\ (1) 

with energy E — E n + h 2 {k 2 x + ky)/2m*, where m* 
stands for the electron effective mass. E n and x«( z ) are 
the eigenvalues and eigenstates of the one-dimensional 
square-well potential problem, with the following prop- 
erties: 

X „(z)~e- Q "- Z (2) 

for z > L/2, assuming the z origin at the center of the 
film, and: 

a n = (V2m*/«) • V^K, (3) 

where the zero of energy is taken at the vacuum level. 
In the three-dimensional problem, E n coincides with the 
energy E of the subband state n at k\\ = 0, measured 
relative to the vacuum level. 

All states ip^k k belonging to subband n, with en- 
ergy E > E n , are thus characterized by the same decay 
length l/a n ~ l/\/—E n , for a given film thickness L. If 
Ef is located between the levels E n and E n+ \, the dom- 
inant decay length of the electronic states at Ef is l/a n , 
namely the decay length of the highest occupied band 
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state at feii = 0. With increasing width L of the well, the 
energies of the quantum-well states decrease with respect 
to Ep) hence, the decay length A of the electronic den- 
sity at Ep first decreases as l/2a„ ~ l/s/—E n , until the 
next quantum- well state at T crosses Ep, at which point 
A increases to the new value l/2a„+i ~ E n+ \. The 

decay length then decreases again with increasing L, dis- 
playing systematic oscillations with L. Considering next 
a discrete number of atomic layers, the periodicity of the 
crossing of the Fermi energy may be derived from the 
Bohr- Sommerf eld rule, which for Mg(0001) yields a peri- 
odicity of 7.7 monolayers (ML) Q. 

From the model description, we expect thus oscillations 
in A, with local maxima occurring when the highest occu- 
pied quantum- well states at Y is closest to Ep. A more 
realistic description of the system, however, is clearly 
needed to confirm the robustness of this behavior, and 
to quantify the effect. In particular, electron interaction, 
atomic orbitals and bonds, surface states, and the pres- 
ence of the substrate, are all factors which are expected 
to influence the decay length A. We have therefore exam- 
ined the effect of quantum confinement on A by means of 
ab initio calculations performed for epitaxial Mg(0001) 
films on W(110) and also, for comparison, for the corre- 
sponding unsupported Mg(0001) films in vacuum. 

First-principles results. — The calculations were per- 
formed within density functional theory, using the 
Pcrdcw-Burkc-Ernzerhof exchange-correlation functional 
0, Troullier-Martins pseudopotentials @], and a plane- 
wave basis set. The Mg pseudopotential was generated 
in the atomic configuration: 3s l s 3p°3d 0/2 , using as core- 
radii cutoff: r S:Pi d = 2.8 a.u.. For W we treated the 6s, 
6p, and 5d orbitals as valence states, using the same pa- 
rameters as in Ref . Q • The overlap between valence and 
core electrons was accounted for using the non-linear-core 
correction to the exchange-correlation potential 0- A 
kinetic-energy cutoff of 49 Ry was used for the plane-wave 
expansion of the electronic orbitals of the Mg/W systems. 
For the isolated Mg(0001) slabs we used a kinetic-energy 
cutoff of 14 Ry. The Brillouin-zone sampling was done 
with a Monkhorst-Pack (MP) grid |9j, and we used a 
Gaussian smearing of the electronic levels of 0.02 Ry to 
determine the Fermi energy. The films were modeled us- 
ing slab geometries in supercells. For the Mg films on W, 
we considered slabs containing 7 W(110) ML, terminated 
by 5 to 12 Mg(0001) ML on one side, and by 2 Mg (0001) 
ML on the other side. The Mg bilayer was introduced in 
order to avoid the presence of an electric field — due to 
the different work functions of W and Mg — in the vac- 
uum regions separating the periodic images of the slab. 
We used vacuum regions with a minimal thickness of 20.1 
A, in order to obtain well-converged values of A. 

Experimentally, above 3 ML, the Mg films are known 
to grow epitaxially on W(110) [j], with lattice parameters 
corresponding to their bulk values fiof . The mismatches 
between the experimental in-plane lattice parameters of 



Mg(0001) (a M g = 3.21 A, 6 Mg = V^Mg) and W(110) 
(a w = 3.16 A, b w = V2a w ) are 1.6% and - 20% 
|4|. To model such epitaxial systems, a commensurate 
interface atomic structure is needed in the calculations. 
The inclusion of dislocations would require prohibitively 
large lateral dimensions of the supercell, and the details 
of the atomic structure at the interface are unknown in 
any case. In order to simulate an unstrained Mg film on 
W(110), we thus elected to laterally strain the W(110) 
slab to the in-plane lattice parameters of Mg(0001). The 
epitaxial alignment was made by positioning the atoms 
of the first W(110) layer, adjacent to the Mg, in the con- 
tinuation of the Mg (0001) hep lattice. The three outer- 
most layers of the Mg film were relaxed, while the other 
interlayer spacings were kept frozen at their bulk posi- 
tion (using the theoretical values of the bulk lattice con- 
stants: OM g = 3.19 A, CM g = 5.18 A, and aw = 3.21 A). 
It should be noted, however that the Mg(0001) surface 
relaxations are small [llj,ll2j, and were found to have a 
negligible impact on A. The self-consistent calculations 
were carried out using a (20,20,1) MP grid. For the local 
density of states, we used a (48,48,1) grid centered at T, 
with a Gaussian level smearing of 0.005 Ry. The decay 
length was derived from a fit, assuming an exponential 
decay of the local density of states at distances beyond 
~ 2.15 A from the outermost atomic plane. 

In Fig.^ w e report the calculated decay length A of the 
Mg films, as a function of film thickness. The results are 
shown for both the Mg films on W and the isolated Mg 
films in vacuum. The behaviors are very similar in the 
two cases. The decay length exhibits a pronounced oscil- 
lation, with a maximum at 9 layers and a minimum at 6-7 
layers. The presence of the tungsten substrate tends to 
reduce the amplitude of the variation of A, from 17% to 
10%, but has no significant impact on the position of the 
extrema. The densities of states at T of the Mg (0001) 
films on W(110) are displayed in Fig. [3 With increasing 
film thickness, an unoccupied quantum- well state crosses 
the Fermi energy at ~ 9 layers. This coincides with the 
maximal decay length A found in Fig. ^ consistent with 
the model prediction. In Fig. |2J we also reported the 
energy positions of the electronic states at V obtained 
from the calculations for the isolated Mg slabs in vac- 
uum. Also in this quantum-well state is found 
to pass through the Fermi energy at ~ 9 layers, which 
corresponds to the largest A obtained in Fig. ^ We note 
that we have also examined the variation of the work- 
function of the Mg films on W with the number of Mg 
layers 0], as quantum size effects can be expected on 
this value ^ij. For the Mg films on W, however, the 
calculated changes in the work function are found to be 
very small, namely a variation of 0.05 eV in the range 
5-12 layers. With increasing thickness, the workfunction 
first decreases, in the range 5-8 layers, from 3.74 to 3.69 
eV, and then increases, in the range 8-10 layers, from 
3.69 to 3.72 eV, and then saturates at 3.72 eV for higher 
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FIG. 1: (Color online) Calculated decay length in vacuum of 
the electronic local density of states at the Fermi energy of 
Mg(0001) films of various thicknesses on W(110) (disks). The 
decay length of the corresponding Mg(0001) films isolated in 
vacuum is also shown for comparison (diamonds). 



coverage (in good agreement with the experimental value 
of 3.66 eV for Mg O). 

Inspection of Fig. [21 reveals a striking correspondence 
between the positions of the peaks in the DOS of the 
Mg films on W and the positions of the levels of the 
isolated Mg films in vacuum. The largest difference be- 
tween the two sets of energies does not exceed 0.2 eV in 
the range [Ep — 8 eV, Ep + 2 eV] . The states indicated by 
"SS" in Fig. [21 originate from the Shockley surface state 
of Mg(0001). In the isolated Mg(0001) slab in vacuum, 
the Shockley states of the two surfaces strongly inter- 
act, giving rise to a pair of split even and odd states — 
with respect to the center-of-slab reflection plane. The 
resulting splitting is as large as 1.4 eV for 5 Mg layers 
and 0.6 eV for 12 Mg layers. In the presence of the W 
substrate, these states persist as strong surface/interface 
resonances, with maxima in the probability density on 
both the outermost (surface) Mg layer and the inner- 
most (interface) Mg layer. The ratio of the probability 
density on the surface relative to the interface layer of 
the high- (respectively, low-) energy state SS tend to in- 
crease (decrease) with increasing film thickness, and is 
~ 1.3 (~ 0.85) in the 12-layer case. In the presence of 
the W substrate, some of the quantum-well states, and 
in particular those in the range [Ep — 6 eV, Ep — 3 eV] 
0,0, remain fully localized within the Mg films — similar 
to the confined quantum- well states obtained from ab ini- 
tio calculations in Ag films on Fe 0] , or sharp Fano res- 
onances, located mostly within the Mg film. Other Mg 
quantum- well states, instead, and in particular the states 
corresponding to the feature which crosses Ep, and is 
indicated by the arrow in Fig. [21 become broader reso- 
nances, corresponding to quantum- well states displaying 
an increased probability density within the W substrate. 
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FIG. 2: (Color online) Calculated density of states at Y 
of the Mg (0001) films on W(110) (solid lines). The film 
thickness increases from 5 to 12 atomic layers (bottom to top 
curve). The densities of states have been convoluted with a 
Gaussian of width 0.2 eV. The symbols (diamonds) show the 
energies of the electronic levels at F of the isolated Mg(0001) 
films in vacuum. The states indicated by "SS" originate from 
the Shockley surface state of Mg(0001) (see text). The zero 
of energy corresponds to the Fermi level. 



The resonant character of the quantum-well states which 
cross the Fermi energy tends to smoothen the jump in A 
observed for thicknesses between 7 and 9 layers. 

Hence, in spite of the presence of the surface/interface 
states, and the resonant character of the quantum-well 
states crossing the Fermi energy, the qualitative behav- 
ior of A, inferred from the particle-in-a-box model, is fully 
supported by the ab initio calculations. The variation in 
the decay length we predict from the ab initio calcula- 
tions is expected to have a direct, exponential impact of 
the electron transfer rate by tunneling — from the metal 
to the O2 molecule — which has been proposed to con- 
trol the initial sticking of the oxygen molecules impinging 
on the surface, via the attractive image charge potential 
on the ionized O^ molecule [3}- Assuming a transfer 
rate by tunneling proportional to e~ d / A , with d the dis- 
tance between the metal surface and the center of mass 
of the molecule, and considering, e.g., a distance d w 3.5 
A — within the expected physisorption range of the O2 
molecules a 10 % variation in A produces a 100 % 

change in the transfer rate. Such a charge is of the order 
of magnitude of the experimental change in the oxidation 
rate at low O2 exposure 0. 

The peak positions of the calculated DOS in the en- 
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ergy range [Ep — 4 eV, Ep], in Fig.|3 compare well with 
the near-normal-incidence PE measurements Q], except 
for a systematic shift, by about +2 Mg ML, of the theo- 
retical spectra with respect to the measured spectra. In 
our calculations, the first quantum-well state with energy 
higher than the surface/interface states SS, enters the 
occupied-state spectrum at a Mg thickness of 8-9 ML. 
In the experiment, this state is found, instead, to en- 
ter the valence-band spectrum at a nominal thickness of 
6-7 ML. The same shift is observed between the calcu- 
lated maximum/minima of A and the experimental max- 
imum/minima of the reactivity. Such a shift is probably 
due to the use of a commensurate, dislocation-free inter- 
face atomic structure, in the calculations, to describe a 
heavily lattice-mismatched epitaxial system. The pres- 
ence of strained Mg layers or dislocations at the interface 
could result in the entrance of a quantum-well state at a 
smaller value of the number of layers of the Mg film. 

The changes in A we predict should be observable, e.g., 
by STM near terrace steps on a Mg surface character- 
ized by regions of different film thicknesses Q. Also, He 
scattering experiments have indicated an apparent step- 
height oscillation in the layer-b y-la ver growth of Pb on 
Ge due to quantum-size effects |l8j . Such changes were 
recently shown to be due mostly to a displacement of the 
topmost layer charge density jl9j , in contrast to previous 
interpretations in terms of a displacement of the surface 
atomic planes. This recent analysis corroborates thus the 
present finding on the behavior of the decay length of the 
Mg films. We note that the changes in the decay length 
of the density of states near the Fermi energy may also be 
related to a recent observation of quantum-size effects on 
the chemisorption properties of Cu(OOl) thin films [20| . 
In addition, layer-KKR (Korringa-Kohn-Rostoker) cal- 
culations |2l| have indicated that the lifetime of negative 
ionic states of molecules, adsorbed on supported metal 
thin films varies with the thickness of the film, through 
coupling to quantum- well states. This effect (not yet 
measured, to our knowledge) was associated, in the cal- 
culations, with oscillations in the amplitude, at the posi- 
tion of the molecule, of the density of empty states above 
Ep. This calculated behavior, consistent with the trend 
we find for A, also supports the interpretation we propose 
for the reactivity changes. 

In conclusion, on the basis of ab initio calculations for 
Mg(0001) ultrathin films on W(110), we have shown that 
the decay length in vacuum of the electronic local density 
of states at the Fermi energy exhibits substantial oscilla- 
tions with the film thickness. The decay length is maxi- 
mal, as a function of film thickness, when a quantum-well 
state passes through the Fermi energy. The changes we 
predict in the decay length are expected to have a ma- 
jor impact on the tunneling rate in the electron transfer 
mechanism, which is believed to control the initial stick- 
ing of O2 on the Mg surface. We therefore propose that 
the experimental tuning of the reactivity is due to quan- 



tum oscillations in the electronic density decay length, 
which should be observable by STM. 
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